Entanglement generation and protection by detuning modulation 
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We introduce a protocol for steady-state entanglement generation and protection based on de- 
tuning modulation in the dissipative interaction between a two-qubit system and a bosonic mode. 
The protocol is a global-addressing scheme which only requires control over the system as a whole. 
We describe a postselection procedure to project the register state onto a subspace of maximally 
entangled states. We also outline how our proposal can be implemented in a circuit-quantum elec- 
trodynamics setup. 

PACS numbers: 03.67.Mn, 03.67.Pp, 85.25.Dq, 42.50.Pq 



I. INTRODUCTION 

The controllable generation of entangled states has 
triggered a considerable amount of interest in the physics 
community. In particular, within cavity-quantum elec- 
trodynamics (cavity-QED) contexts, the achievement of 
two-qubit entanglement has seen a flourishing of propos- 
als. It has been suggested to set entanglement between 
two remote atomic qubits by using the cancellation of 
which-path information relative to spontaneously emit- 
ted photons The resonant interaction of a cavity 
field with two qubits has revealed a striking entangling 
power even when the field is prepared in an incoherent 
state [2]. In a similar setup, entanglement can be created 
through the continuous detection of the field leaking from 
the cavity containing the atoms 0. However, qubit en- 
tanglement can also be set in a regime of dissipative dy- 
namics, where the system at hand interacts with a struc- 
tured environment [3] and the conditions for entangle- 
ment generation between subsystems undergoing purely 
dissipative dynamics have been studied 0, 01 ■ Strategies 
for the effective engineering and simulation of such envi- 
ronments have subsequently been envisaged [7] . Most re- 
cently, schemes have been proposed for efficiently induc- 
ing discrete-variable entanglement in a bipartite system 
by transferring the correlation properties of a continuous- 
variable state Q. 

It is easy to recognize the importance of protocols that 
are able to reliably protect correlations, once they are es- 
tablished, from the unavoidable spoiling effects of deco- 
herence and decay. This has led to proposals of passive 
as well as active schemes [l]. The use of decoherence- 
free subspaces is the prototypical example of passive 
strategies, where the correlations set in a system can be 
protected by choosing a proper encoding of the infor- 
mation [l^, [llj. More recently, it has been suggested 
to use macroscopic quantum jumps as a tool to cre- 
ate an entan gled state of two qubits, preparing it in 
a dark state 'l^. One could also consider generalized 
dynamical-decoupling (bang-bang) schemes [13] to ac- 
tively cancel the effect of the environment on the system 
of interest. These strategies are appealing and intrigu- 



ing from a theoretical point of view and proof-of-principle 
experiments have been performed, for example in a solid- 
state setup 114]. However, they are still far from complete 
and too demanding for the current state of the art. 

In this paper we find a strategy, based on dissipative 
qubit-bus dynamics, enabling the simultaneous genera- 
tion of two-qubit entanglement and protection against 
both qubit and bus losses. No structured-bath engineer- 
ing is required in our scheme. The protection is achieved 
by using a simple global addressing of the register, a fea- 
ture which relaxes the usually assumed requirement of 
single-qubit addressing at the center of many dynamical- 
decoupling protocols and brings our scheme closer to 
experimental feasibility. While a quantitative analysis 
is deferred until later, here we briefly provide an intu- 
itive picture of the mechanism behind our proposal. It is 
known from the study of the so-called Dickc model fT5| 
that a bipartite qubit system, prepared with the qubits 
in their excited state and exposed to the fluctuations of 
a common reservoir, soon decays via the channel given 
by the symmetric state \s) = (l/\/2)(l01) -f JlO)) (with 
]0) and ]1) the single-qubit logical levels) into the total 
ground state of the qubits. Thus, there is a transient 
period in the dynamics of the two qubits when a max- 
imally entangled component is involved in the state of 
the system. However, the system does not exhibit en- 
tanglement because of a competition between the fad- 
ing symmetric state and the increasingly populated col- 
lective ground state. In order to reverse this situation, 
the influence of \s) has to be emphasized and stabilized 
with time. The protection from environmental effects is 
then achieved by this stabilization and also by simulta- 
neously inducing a relative phase between the qubits. In 
proper conditions, this results in part of the population of 
the symmetric state being moved into the antisymmet- 
ric state \a) — (l/\/2)(]01) — JlO)), which is decoupled 
from the decay mechanism (it is a subradiant state) [15| . 
In this paper we show that an entangled steady state is 
produced by modulating a detuning between the com- 
mon bus and the register. The results depend on the 
temporal profile of the modulation, which represents a 
dial that can be turned so as to span a specific sector 
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of entangled steady states. Finally, we address an ac- 
tive protocol based on postselection, which projects the 
state of the qubits almost completely onto the symmetric 
state, thereby achieving nearly maximal entanglement. 

It is important to stress the differences between a 
bang-bang scheme [3] and our own one. Although both 
scheme ultimately rely on the control of the interaction 
between the register and the environment, the two ap- 
proaches are intrinsically different. Bang-bang schemes 
keep the state of interest unchanged throughout the evo- 
lution by effectively decoupling it from the environment. 
Our protocol however produces entanglement through 
the exploitation of purely dissipative dynamics. The pro- 
tection of the entanglement from the influences of the 
external world is provided by the development of a sub- 
radiant behavior of the system due to the detuning mod- 
ulation. Moreover, in a bang-bang protocol the timing is 
set by the fast switching rate of a control field, which is 
given by the inverse of the coupling strength between the 
system and the environment. Our scheme, on the other 
hand, is based on a weak-coupling regime between the 
register and the bus, which sets a slower time-scale than 
bang-bang schemes. This reflects not just a technical di- 
versity, but is a manifestation of two almost complemen- 
tary ways of designing protection from an environment. 

The remainder of the paper is organized as follows. In 
Section [IT] we address the Bloch equations derived from 
the qubits' reduced master equation, in a weak-coupling 
regime and first Born-Markov approximation. This de- 
scribes the dynamics of the qubits interacting with a 
leaky bus mode exposed to a bosonic multi-mode reser- 
voir at thermal equilibrium with temperature T. The 
Bloch equations (as well as the master equation) fully 
account for a time-dependent detuning, which modulates 
the couphng between the qubits and the bus. Section Hill 
explains in some detail how single-qubit addressing is not 
required in our scheme. This paves the way toward set- 
ting our proposal within the context of active (but sim- 
ple) global-addressing scenarios. In Section [IV] we de- 
scribe the postselection protocol to improve the amount 
of entanglement set between the qubits by a nearly per- 
fect projection of their joint state onto the symmetric 
state \s). Finally, Section |V] describes a physical system 
to embody our proposed scheme. We show in some detail 
how a circuit-QED setup of two superconducting charge 
qubits incorporated in a microwave planar stripline res- 
onator can be used. This offers practical advantages com- 
pared to the standard cavity-QED setup. The details re- 
garding the derivation of the qubits' master equation are 
given in the Appendix. 



II. THE SYSTEM AND DETUNING 
MODULATION PROTOCOL 

We consider two qubits, labelled 1 and 2, interacting 
with a one-sided single-mode cavity field described by 
the bosonic creation (annihilation) operator (d) . Each 
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FIG. 1: (Color online) Two qubits interact with the same cav- 
ity mode (in contact with a thermal reservoir, decay rate k) 
acting as a common bus. An inhomogeneous global potential, 
setting a time-dependent detuning is also shown. 



qubit is characterized by a transition frequency E'^ and 
interacts with the cavity field (frequency lOq) via a dipole 
interaction with strength g. We assume a leaky cavity ex- 
posed to a multi-mode bosonic environment. The cavity 
mode is in a thermal state with temperature T and aver- 
age photon number h = (e^'^" — where f3~^ = KbT 
and Kb is the Boltzmann constant (we use units such 
that h= 1). The strength of the field mode-external bath 
coupling is given by the cavity decay-rate k. Within the 
rotating wave approximation, we model each qubit-field 
interaction as Hcj ~ g{a^aj -t- h.c.) {j = 1,2), where 
aj = ((Tj^)^ = |0)j-(l|. The total energy of the system is 



Hsys = uJc a^a + ^ + H,j] (1) 

with (t| the z— Pauli matrix of qubit j. We have con- 
sidered an implicit time dependence of the single-qubit 
transition frequency. The detuning between the j—th 
qubit and the cavity field mode is indicated as Aj(t) — 
Ej{t) — ujc- By introducing the transverse-mode energy 
decay rate 7 and assuming k <C Wc, the evolution of the 
qubits-bus system is described by the master equation 
Eq. ([5]) given in the Appendix. Our system is sketched 
in Fig. m 

As we are interested in just the dynamics of the 
qubits, we adiabatically eliminate the bus mode and de- 
rive the corresponding Bloch equations. This is straight- 
forwardly done by deriving a reduced master equation 
for the qubits only and projecting it onto the basis 
{|f= 11) , |s) , \a) , 11= 00)}i2. We refer to the Appendix 
for the details of the adiabatic elimination. Here we con- 
centrate on the form of the Bloch equations relevant to 
our work, where we assume the initial state |T)i2 is pre- 
pared and, for the sake of simplicity, we consider the case 
of a detuning in amplitude smaller than the cavity decay 
rate. Using the notation gij = i2(i|p|i)i2 (*,i =1,s,a,l) 
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FIG. 2: (Color online) Fidelity of the two-qubit state g for the 
full-resonant case with the symmetric state |s)-^2 (solid line) 
and the antisymmetric one |a)]^2 (identical to zero) against 
the dimensionless interaction time r : 



: g't/K. 



(2) 



and G1,{kn) — qj + ^ they read 

5t£»TT = -4G}(n)£<n + 2GSJ(n)(e,, + gaa), 
dtQss = -Gi{2n)[cosA{t)gss - ismA{t)gsa + h.c' 

+ 2G}(n)(e„ - Qss) - 2G°(n)(^^,, - 

+ 2cos A(t)[G°(n)gn + G°(n)f?xJ, 
dtQsa = -2G\{2n)gsa + ismA{t)[2G'i{n)g^^ 

+ 2G°{n)gu - G'i{2n)igss + Qaa)]. 

The equation for gaa is given by that for gss with gss — > 
gaa and A(i) ^ tt — A(t). The equation for gas is the 
hermitian conjugate of the one for gsa and the equation 
for gii is found from g^ = 1 — — gss — gaa- We have 
taken Ai{t) > 0, with A2{t) = and Ai{t) = A{t) for 
ease of notation. In Section UTTl we consider the general- 
ization of this situation to two-qubit detuning and in the 
Apppendix we provide the form of the reduced master 
equation. Moreover, we stress that the absence of terms 
like g^a, g^s (and analogous) is due to the specific choice 
of the initial state. In particular, if any coherence is ini- 
tially present in the qubit state, the set of Eqs. ^ must 
be complemented by a second closed system of Bloch 
equations which can easily be derived. 

From Eqs. ([2]) the initial state \1)i2 evolves into 



C{g), we relax the max condition. Thus, in the follow- 
ing plots, entanglement is present only when C{p) > 0. 
Moreover, we address a physically interesting situation 
by considering the initial state |T)i2 which, as our sys- 
tem is formally equivalent to a Dickc model decays 
toward the ground state ||)-^2 on a time-scale which is 
faster than the single-qubit relaxation time. Thus, for 
the steady state of the system, no entanglement is ex- 
pected between 1 and 2. As we stress later, this choice 
for an initial state is dictated by the global addressing 
context of this work. The state |T)i2 can be prepared via 
a global potential and without local control [17| . 

As time passes, the superradiant state is rotated to- 
ward a mixed state. While still exhibiting no quantum 
correlations (it never violates the Peres-Horodecki sepa- 
rability criterion [3), the state nevertheless has a good 
projection onto |s)-^2 ^'^d no contribution from the an- 
tisymmetric state \a)-^^2■ This is shown in Fig. [5] where 
the fidelities J-j = i2{j\g\j)i2 = Qjj {j = s, a) are plotted 
against the dimensionless interaction time t = g^t/ k. At 
T ~ 2.5, gss ^ 0.37 (with gaa = 0), though the state 
is still mixed and separable due to the non-zero value of 
the fully polarized states |t)i2 and \i)i2- This suggests a 
minimization of the infiuence of these components on g. 
The idea behind our proposal is to exploit this fact and 
change the dynamical evolution of the qubits, by intro- 
ducing a detuning A, so as to induce an evolution with 
an initial state which is no more |T)i2 but p(r — 2.5). 
Thus, at T = 2.5 an external potential is switched on 
and changes the transition frequency Ei(t). The dynam- 
ics of the qubits are therefore described by Eqs. ^ (no 
component outside the subspace encompassed by Eqs. ^ 
is present in the new initial state) and the concurrence is 
plotted against r in Fig. [3] (solid line), where A = lOg^/n 
with g — 0.3k, so that the adiabatic condition is fully re- 
spected. We find C > 0.3, stable at the steady state of the 
qubits. This interesting result can be compared with 
which does not achieve a stable entanglement. The plot 
results from a transition between the entanglement func- 
tion of the full resonant case before t = 2.5 and that 
of the single-qubit detuning described above for r > 2.5. 



z2 y)i2{j\ + iSas\a) 12(^1+ h.c). (3) 

j=1,a,s,l 

In order to illustrate the basic features of our proposal, 
we consider 7 = fi = for the moment. These parame- 
ters will be re-introduced later on. As a measure of the 
entanglement in the bipartite mixed state of qubits 1 and 
2, we use the concurrence (l6j which, for a bipartite state, 
can be calculated as C{g) = max[0, -^0! — J2i V^] with 
g = £'('X'^^icr0£'*(<Xi^^iO'^) (g* is the complex conjugate 
of g) and ai > ai{i ~ 2,3,4) are the eigenvalues of g. 
For a qubit state as Eq. ([3]), we have 



C{g) = max[0,2(|eio,oi| - V^nQll)] 



(4) 



with 2£>io,oi = Qss + Qsa - Qas ^ Qaa- In Order to gain 
as much information as we can about the behavior of 




FIG. 3: (Color online) C(p) against r for the initial state |T)i2 
and A(r) = lO0(r — ro). We have considered ro = 1.5, 2.5, 3.5 
(dashed, solid and dotted line respectively). The inset shows 
the behavior of the detuning functions. 
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This situation is equivalent to taking A(t) = 108(t— 2.5) 
(in units of / k), where 6(t — tq) is the Heaviside func- 
tion. The steady-state value of the entanglement turns 
out to dependent weakly on the amplitude of A(t) but 
strongly on tq. For instance, in Fig. [3] we show the re- 
sults of small deviations from the case considered above 
by plotting the concurrence relative to tq = 1.5 (dashed 
line) and tq = 3.5 (dotted line), which give rise to smaller 
steady-state entanglement. This results from a smaller 
Qss{tq) component in the two-qubit state and a disad- 
vantageous competition between \s)^2 1*^)12' which 
lowers the entanglement. This is strikingly exemplified 
by increasing tq by one order of magnitude. In this case, 
the switching on of the detuning occurs in correspondence 
to a state of the two qubits which is mainly decayed to 
Ii>i2 (^'u(To = 25) = 0.999501). 

It is worth stressing that, even though our choice for 
g/ K may seem to put the above example at the boundary 
of the applicability of the adiabatic elimination, we have 
checked that for g/ k, ^ 0.1, no significant change occurs 
in the entanglement generation process [30l |. In principle, 
the true evolution of the system, obtained by numerically 
solving the complete master equation (without adiabatic 
elimination) should be compared to the situation here at 
hand. However, this is in general a very hard task, which 
goes beyond the scopes of the present work. Nevertheless, 
by checking the effects of values for g / k which are largely 
within the validity of the adiabatic elimination, we can be 
sure of the validity of the above approach. The system 
is flexible enough to tolerate a less strict ratio of the 
different time-scales involved in the problem. 

The appearance of steady-state entanglement can be 
shown clearly by the behavior of the density matrix ele- 
ments. A careful analysis reveals that the qubit entan- 
glement is due to the presence of \s)^2 ^^"^ I^)i2- '^^^ 
calculation of the fidelities Ts,a in presence of the detun- 
ing modulation, reveals that as soon as the detuning is 
switched on, an \a)^2 component is developed, which af- 
ter a transient period, stabilizes to a steady state value. 
This stabilization holds also for the \s)i2 component, 
with Qgg ^ Qaa- Thcse behaviors can be seen in Fig.|3]for 
the Heaviside function with tq = 2.5. However never 
vanishes, thus affecting the entanglement. This clarifies 
the mechanism behind the entanglement generation and 
protection. Without a detuning modulation, the system 
would never develop any subradiant behavior {i.e. the 
overlap with |a).|^2 always be zero). This is not true 
for a modulated situation, where the dynamical condi- 
tions are changed. Once the system has decayed into 
an incoherent superposition of |s)i2' It) 12 ^^'^ Ii)i2 (^^^^ 
r < To), one qubit acquires a relative phase with respect 
to the other one due to the detuning, which results in 
the development of a subradiant component. The steady 
state entanglement is the result of a competition between 
||)-^2! the antisymmetric and the symmetric component. 

The assumption of a Heaviside function regulating 
A(t) is not critical. In order to relax this assump- 
tion, we have checked the results corresponding to a 
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FIG. 4: (Color online) Fidelities Ts,a for a detuning modu- 
lation strategy based on a Heaviside function with tq = 2.5. 
The solid line is for Ts, the dotted one for Pa- At r = to, 
an antisymmetric state component is developed. For large 
r, both the antisymmetric and symmetric state fidelities are 
stabilized. 



smooth raising edge of the detunings given by the func- 
tion A[l + e^'''^'^"~'^']~2 with A an amplitude. For proper 
choices of A and 6, this is a slowly rising function (with 
respect to the time-scale set by k, see the Appendix) 
producing a concurrence which differs from the result ob- 
tained for a Heaviside function by less than 1%, as shown 
in Fig. El for 6 = 3, tq = 2.5 and A = 10. Thus, for defi- 
niteness and in order to simplify the calculations, we as- 
sume a Heaviside profile of the detuning rising edge [sH . 

It should now be clear that a detuning function with 
a single rising edge represents the best choice. This can 
be confirmed by considering the value of C as a function 
of the time width of a single square pulse. As already 
stressed, the turning on of the detuning corresponds to 
the maximization of the symmetric state component and 
the introduction of an \a)^2 component in g. If the detun- 
ing is switched off after a time 5r, the symmetric com- 
ponent quickly goes to zero (together with any correla- 
tion betwen \s) and \a)) while the subradiant part is pre- 
served. This achieves a non-zero stationary entanglement 
which nevertheless is reduced with respect to the case of 
a Heaviside function. Indeed, in the above conditions, at 




FIG. 5: (Color online) Comparison between the concurrence 
obtained by using a Heaviside function lO0(r — 2.5) (dashed 
line) and the smooth detuning function 10[1 -|- e"^^'^^"^']^ 2 
(solid line) . The inset shows the time behavior of the detuning 
functions. 
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the generic time Tg , the steady-state entanglement quan- 
titatively corresponds to the fraction of the antisymmet- 
ric state being present in p, as can be immediately seen 
by considering a state like g — X]j=t s a [-^i b)i20l- ^"-"^ 
Te — tq 3> 5t, the population of |t)]^2 ^'^'^ 
the symmetric state fraction. That is A| ^ = so that 
C{q) = Aa- The entanglement is stable, even though 
small, due to the sole presence of the subradiant compo- 
nent, developed at the switching-on of the pulse. As soon 
as St becomes larger than Te, making the symmetric state 
fraction (and its correlations with |a)j^2) non- negligible, 
the entanglement is not only stable but also reaches the 
asymptotic value corresponding to Fig. [31 as can easily be 
seen. As Te is increased, this behavior holds for a larger 
St, demonstating that a large steady-state entanglement 
is achieved only for a step-like function (as mentioned be- 
fore, the raising edge functional behavior is irrelevant). 

A further case can be considered, namely a peri- 
odic modulation. However, this modulation implies the 
switching on/off of the detuning at instants of time that 
correspond to smaller fidelities of the state g with the 
symmetric state. For instance, in Fig. [S] we consider the 

case of a square wave Ast„(r) = 10 X]n=i(~l)"^^0(''' ~ 
2.5n), which produces N/2 square pulses of amplitude 10 
(units of g^/zi). No entanglement results from this detun- 
ing modulation strategy, as the switching of the detuning 
introduces a jagged-drop in the fidelity compared to the 
situation in Fig [H While the on part of the square wave 
always corresponds to a slow increase of C, the off part 
results in a larger decrease, resulting in an overall pull 
down of entanglement. 

So far, only the case with no spontaneous emission and 
a zero-temperature bath has been considered. In order to 
include the effects of 7,n ^ 0, we need to solve Eqs. ^ 
for 7, n ^ and relate them more closely to a physical 
setup that will implement our protocol. Although more 
detail will be given later, here we mention that the situa- 
tion considered is such that 7 ^ c/^/k and n 1, which 
are realistic conditions in several physical systems such 
as circuit-quantum electrodynamics of superconducting 
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FIG. 6: (Color online) Concurrence (red curve, left vertical 
scale) superimposed on the detuning function Asm(r) (black 
curve, right scale) against r for the qubits initially prepared in 
the superradiant state |T)i2- In this plot we have considered 
= 16. No entanglement is ever present between the qubits. 



charge qubits integrated in microwave cavities [19|, |20|. 
We will postpone discussion of the order of magnitude 
of these physical parameters until Section |Vl To fix the 
ideas and to be as close as possible to physical reality, 
we consider n = 0.06 and 7 = 10~'^k. Moreover, in tack- 
ling this analysis, we find it convenient to refer to the 
computational basis {|t) , |10) , |01) , ||)}i2. The essen- 
tial features of the previously considered case still hold, 
with the Ts function also maximized at t = 2.5. The 
main effect of the non-zero thermal photon number is a 
reduction in the steady-state entanglement value. This 
is due to n ^ 0, which reduces the coherence pio.oii thus 
preventing the state from mimicking the symmetric state. 
On the other hand, 7 7^ introduces a second time-scale 
in the system, which results in a slow decay of both the 
populations of states |10)-^2 1 |01)i2 ^^'^ '^^ coherence 
Pio.oi- This accounts for an overall decay of C, as shown 
in Fig. [71 Despite the decrease found for non-zero 7, the 
entanglement still remains as large as 0.2 for interaction 
times up to t = 100. This corresponds to an effective 
entanglement protection from both thermal behavior of 
the bosonic bath and the qubits' spontaneous decay. 

III. DOUBLE DETUNING AND GLOBAL 
ADDRESSING 

In this Section we address the question of whether 
the assumption of a single-qubit detuning modulation is 
critical to the proposed protocol. The answer to this 
provides an effective justification and an a posteriori 
motivation for the analysis conducted so far. In or- 
der to do this, we re-consider the physical system de- 
scribed in Section [TTl (see the Appendix) with the in- 
clusion of the second qubit detuning A2 (t) . This modi- 
fies the dynamics of the system in such a way that the 
Bloch equations ^ remain unaltered with the replace- 
ment A(t) A(t) = Ai{t) - A2(t). Thus, considering 
both the qubits as detuned in time is equivalent to just 




FIG. 7: (Color online) Concurrence against r for a modulated 
detuning with g — 0.3k, 7 = 10~^«; and n — 0.06. The time 
axis has been extended in order to show that a concurrence 
close to 0.2 is present for r up to 100. The slow entanglement 
decay is due to the presence of a non-zero 7, while the smaller 
steady-state value is largely due to the thermal nature of the 
bath (n / 0). 
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considering the energy of qubit 1 being modulated with 
an effective time-dependent detuning A{t). That is, the 
analysis conducted so far is perfectly general and there is 
no limitation in considering a single-qubit modulation as 
this case encompasses rigorously the most general situa- 
tion of double detuning. Obviously, the Aj{t)'s cannot 
be chosen arbitrarily, in general. 

This result has two main implications. The first is that 
in order for the protocol we have described to be effective, 
we must consider detuning functions which are opposite 
in sign {i.e. one detuning has to be positive, the other 
negative). The second point is pragmatically relevant as 
we can now put our scheme within the context of global 
addressing protocols [2l|. Indeed, it should be clear after 
the above discussion that the realization of the detuning- 
modulation protocol simply requires the appropriate set- 
ting of a potential which addresses both the qubits in the 
correct way (increasing the energy spacing of qubit 1 with 
respect to the resonant value coc and reducing the spac- 
ing of qubit 2). No single-qubit addressability is needed, 
which considerably reduces the experimental efforts for 
the implementation of the scheme. It is not necessary to 
require a strongly focused potential applied to just one of 
the two qubits and having no effect on the dynamics of 
the other one. In order to fix this idea, one can consider 
a global magnetic field inducing a Zeeman-like effect on 
the qubits' energy levels, the shifting being different from 
qubit to qubit because of a gradient in the magnetic po- 
tential (see Fig.[T|). In Section|V]we address the physical 
mechanism responsible for such a shift by considering a 
specific experimental setup that can be used in order to 
implement our proposal. 



IV. ENTANGLEMENT IMPROVEMENT BY 
POSTSELECTION 

As the entanglement is set in a (quasi) steady-state, 
the required degree of control over the system is re- 
duced. With the exception of the choice for the opti- 
mal value of To at which to switch on the detuning, no 
fine time control is necessary in order to properly drive 
the dynamics of the system. Nevertheless, it might be 
desirable in many situations to raise C up to a maxi- 
mal entanglement of one ebit. We have seen an intrin- 
sic limitation in the amount of establishable entangle- 
ment due to the unavoidable presence of the spurious 
population of |4)-|^2- ^ procedure which allows us to 
cut away the unwanted contribution from \i)i2 is rep- 
resented by the postselection of the two-qubit state af- 
ter some detection event. Explicitly, consider the fad- 
ing influence of the component. In the specific case 
here at hand, each time the state of the two qubits is 
not found to be \l)i27 the overlap with |s)-^2 increases, 
improving the entanglement between the qubits. Thus, 
by using the positive-operator-valued-measure (POVM) 
{Ho = Ii)i2(ihni = 1 - Ho} with 1 the identity op- 
erator, we can postselect the state resulting from the 
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FIG. 8: (Color online) Comparison between the entangle- 
ment in the postselected state with and without the detuning- 
modulation protocol (solid and dashed curves respectively). A 
Heaviside modulation is considered in the solid line case. 

qubits not being found in the global ground state. This 
changes g into Qp = Affligfli = N(q - t^ii U)i2(il) 
with M a normalization factor. As gf^f — > 0, this ef- 
fectively results in a projection of the two-qubit state 
onto the subspace spanned by \s,a)-^2 with asymptoti- 
cally i2(s|gp|s)i2 > 0.9. After the analysis in Section 
II, we know that a large fraction of |s)-^2 implies a large 
degree of entanglement. This is witnessed by the entan- 
glement properties of the resulting state, which is repre- 
sented in Fig.lHl The plot represents the amount of entan- 
glement in the postselected state when the measurement 
is performed at the instant r in the evolution of the two- 
qubit state. Both the detuning-modulated (solid line) 
and the full-resonant case (dashed line) are shown. In 
both the cases there is an improvement in the amount of 
stationary entanglement. While the second case has con- 
currence stabilized around 0.4, an entanglement larger 
than 0.9 is obtained in the modulated case. However, a 
full ebit is not possible as the spontaneous emission and 
the thermal effects of the bath spoil the correlations be- 
tween the qubits. Indeed, we have checked that almost 
a complete ebit is achievable if 7 = fi = is considered, 
where the detuning-modulated condition is better than 
the unmodulated the concurrence value always 

remains above the entanglement curve of the unmodu- 
lated case. 



V. PHYSICAL SETUP 

The physical system we suggest to use in order 
to implement our proposal is given by two supercon- 
ducting charge qubits, embodied by Superconducting- 
Quantum-Inferference-Devices (SQUIDs) [12], nano- 
lithographically implanted in a quasi-unidimensional mi- 
crowave stripline resonator [2^. This system offers ad- 
vantages in many respects. First, the qubits are station- 
ary within the cavity, so that the requirement for a fine 
tuning of the transit-time through the cavity (typical of 
microwave cavity-QED implementations) is no more an 
issue. Second, the coupling between the qubits in the 
register and the cavity bus can easily be arranged so as 
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to satisfy the weak coupling regime required by our pro- 
posal. Finally, the manipulability of charge qubits em- 
bodied by SQUIDs allows for a detuning modulation in 
the global addressing fashion depicted in this paper. 

In detail, we assume the charging regime and the low- 
temperature limit [l^l and set each SQUID to work at the 
charge degeneracy point, where the qubits are encoded 
in equally-weighted superpositions of states, having zero 
and one excess Cooper-pair on the SQUID island, namely 
|±)^. = (l/\/2)(|0) ± \2e))j (2e being the charge of a 
Cooper pair). The degeneracy point is set by biasing each 
SQUID with a dc electric field connected to the super- 
conducting devices via the ground plate of the resonator. 
The free Hamiltonian of a single SQUID is thus given by 
(l/2)i;°(0)CT|, with E^j{(j)) the Josephson energy (tuned 
via an external magnetic flux </> piercing the SQUIDs). 
By modulating the magnetic flux, we can change the en- 
ergy separation between the qubit levels, thus setting the 
detunings with respect to the cavity mode frequency. A 
gradient can be incorporated into the external magnetic 
flux so as to realize a configuration of equal and opposite 
detunings for the qubits in our register. 

The microstrip resonator can be modelled as a dis- 
tributed LC oscillator, where C is the capacitance be- 
tween the plates of the stripline and L is the overall in- 
ductance of the device (depending on the length of the 
resonator, typically in the range of 1 cm). In this setup, 
7/k ~ 10"'^ and a cavity quality factor of ~ 100 are con- 
servative assumptions. At Wc/Stt = 6 GHz and T ~ 170 
mK we have n ~ 0.06. The coupling between qubits and 
cavity mode is capacitive and mediated by the electric 
part of the cavity field [l^, [2^. In a second quantiza- 
tion picture, the interaction Hamiltonian can be cast in 
the form of a Jaynes-Cummings model so that Hgys in 
Eq. ([1]) can naturally be embodied by the present setup. 
A Liouvillian description of SQUID-cavity open systems 
has been proven to be rigorous up to temperatures well 
above those assumed in this work. Indeed, for two qubits 
in a stripline resonator, the optical master equation ^ 
can be derived from the Bloch-Redfield formalism, when 
the secular approximation is relaxed and a large number 



of elements of the Redfield tensor are considered 

Two SQUID qubits (size ~ /^m) can easily be accom- 
modated in the cavity far enough away from each other 
to achieve negligible cross-talk (in principle due to direct 
capacitive and inductive coupling). Lithographic tech- 
niques allow us to control, within a few percent, the ge- 
ometric characteristics and the resulting parameters of 
the device. The two qubits can therefore be manipulated 
both simultaneously or independently with two separate 
coils. Due to charged impurities in the vicinity of the de- 
vices, separate calibration at the degeneracy points would 
be required for each qubit. This may be achieved with 
several adjustments to the design of the setup. For in- 
stance, by splitting the ground plate and attaching a gate 
to each part p^ . 

Let us now turn briefly to the description of possible 
ways of implementing the conditional detection scheme 



described in Section IIVI In principle, a measurement 
of the qubits' state can be performed by setting a large 
qubit-cavity field detuning, attaching a detector at the 
output capacitive gap of the stripline resonator and mea- 
suring the shifts induced in the resonance spectrum of a 
probe beam sent into the cavity through the input ca- 
pacitive gap. The dispersive nature of the qubit-cavity 
coupling, which changes the refractive index of the cav- 
ity field mode, determines qubit-state dependent shifts 
in the resonance peak of the probe beam. This allows for 
the non-demolition detection of the qubit state, follow- 
ing the strategy depicted in [23|. However, in order for 
these shifts to be detectable, the change of the refractive 
index has to be larger than the cavity linewidth, a con- 
dition which is hard to match if the bad cavity regime is 
invoked. However, a second strategy is possible, which is 
more suitable for conditions of large detunings between 
the cavity and register and a large cavity decay rate. 
This involves driving a cavity field mode with a coherent 
state \a) (a G C). In the situation of a large qubit-cavity 
detuning, the dispersive dynamics the system undergoes 
is such that the globally unexcited state \Vi\2 becomes 



correlated with the field state |ae^®) [20|, |25l. That is, 
in phase space the coherent state acquires an additional 
phase dependent in general on the ratio On the 

other hand, the symmetric and antisymmetric compo- 
nent of the density matrix leave the coherent state un- 
changed [2^. A homodyne measurement of the cavity 
field provides a distinction between the states of the reg- 
ister and therefore the implementation of the POVM we 
have described. 

As an additional remark, we stress that in this setup, 
at the charge degeneracy point, decoherence due to low- 
frequency modes vanishes at the first order. This allows 
the minimization of the effects of noise sources repre- 
sented by switching charged impurities in the proxim- 
ity of the SQUIDs' islands, which constitute a system of 
bistable fluctuators giving rise to 1// noise [2^. Finally, 
it is worth stressing that due to the qubit-resonator in- 
teraction, the energy levels of our qubits are much less 
sensitive to these charge fluctuations than isolated qubits 
at the optimal working point [l9| . This allows us to ne- 
glect any resulting dephasing effects. 

Finally, it is worth mentioning that the example con- 
sidered in this Section is just one of the possible physical 
setups where our proposal could be implemented. In- 
deed, the formalism we have used in order to describe 
the main features of our protocol is general enough to be 
adapted to various situations. For instance, the case of 
two trapped ions, in the Lamb-Dicke regime and placed 
inside an optical cavity can be taken in consideration [2^ . 
The extension of our analysis to the case of multi-level 
systems composing the register, on the other hand, will 
pave the way to the use of two closely-spaced ensem- 
bles of cold two-level atoms (confined in vapor cells or 
magneto-optical traps). The free-space interaction of a 
laser with the ensembles, each treated as an effective 
A'^/2-spin (where N is the number of atoms in each ensem- 
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ble) and within the rotating wave approximation, pro- 
vides an interaction Hamiltonian which is the generahza- 
tion of our model to iV + 1 systems [29;] . 
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APPENDIX: MASTER EQUATION 

Single-qubit detuning 

In this Appendix we briefly describe the steps taken 
in order to derive the Bloch equations ([2]). The last part 
of the Appendix is dedicated to the case of a double de- 
tuning, as studied in Section [TTTl In an open- system per- 
spective, the dynamics of the cavity field-qubit system is 
described by the master equation 

2 

dtPcl2 = -i[Hsys,Pcl2] + C1[pcl2] +^£j[Pcl2], (5) 

J = l 

where Hsys is defined by Eq. ([T]). Here, we have intro- 
duced the Liouvillian superoperators describing the cav- 
ity decay £c[Pci2] = n{n + l)[2apci2a) - {a^a,pci2}] + 
Kfi^a) pci20. — {aa\ pci2Y\ and the single-qubit sponta- 
neous emission >Cj[pci2], which has the same structure as 

^c[Pci2] for n = 0, K ^ 7 and a ^ &J . 

As we are interested in the qubits' evolution, we now 
adiabatically eliminate the degrees of freedom of the field 
mode. This can be done using a procedure valid in the 
weak-coupling regime g <ti k. In this case, the dynamics 
of the mode interacting with the bath is much faster than 
its interaction with the qubits. The qubits see the cavity 
mode in a steady state pss not affected by the qubit- 
mode dynamics and determined by the statistical prop- 
erties of the environmental bath. In the case at hand, 
Pss is a thermal state with average photon number h. 
By using, standard techniques of quantum optics based 
on second-order perturbation theory and the first Born- 
Markov approximation |27| (see also ref. [Bj), the qubit 
master equation in the interaction picture reads 

2 

dtg = Y^icf^'^ [g] + C][g]) + Ciy^-ig] (6) 



with 

2 

Cf/'^^ [g] = + 1) [2ajga+ - {a+aj, g}] 

(7) 

+ n[2ij+gaj -{&-&+, g}] , 
Cl^^ig] ^{n + ma,_ga+ - {a+a^,g}] 
At) 

H n\2&2 Q'^i ^ {'^i '^2 ^ q}] + h.c, 

K 

and g{t) — 5e~5^(*). Here, we have considered the case 
of E^{t) > coc with E^{t) = E^=LUc (so that A2{t) = 
always) and Ai{t) = A{t) for ease of notation. Eq. ([7]) 
is the Liouvillian describing the single-qubit decay rate 
induced by the coupling to the external thermal bath me- 
diated by the bus. In general, the cavity-induced decay 
rate of atom j depends on the corresponding detuning. 
Within the range of validity our approximations, how- 
ever, the bare decay rates can be used. It is worth stress- 
ing that as the adiabatic elimination procedure does not 
affect the qubits degrees of freedom, the form of ^'jig] 
is left unchanged. However, the adiabatic elimination of 
the cavity gives rise to an effective qubit-qubit interaction 
term. The structures of Eqs. ([7]) and ^ result from the 
assumption that the detuning A(t) is modulated within 
a time-scale slower than the one set by the cavity decay 
rate n. 



Two-qubit detuning 

The introduction of the second-qubit detuning pro- 
ceeds as described in Section |llT] and by re-considering 
Eq. By revising the adiabatic elimination proce- 

dure, the form of the single-qubit decaying terms in 
Eq. ([7]) remain unaltered, while the inter-qubit correla- 

tion term is changed into [g] with g{t) — ge^i^'-*\ 

A{t) = Ai(i) - A2(i) and Aj{t) the detuning of qubit j. 
Evidently, the form of the master equation is invariant 
after the introduction of the two-qubit detuning and the 
dynamics depend only on the relative detuning between 
the qubits. This is due to the structure of Eq. ([5]) in the 

adiabatically eliminated form. In [g], the presence 

of a lowering ladder operator {i.e. aj) is accompanied 
by the rising ladder operator of the other qubit. As a 
lowering operator is associated with the time-dependent 
term e~'^i'^*\ the accompanying rising operator will in- 
troduce so that they always combine to give the 
relative detuning A{t). 
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